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Abstract. In this paper, we prove the existence of two-dimensional, traveling, capillary- 
gravity, water waves with compactly supported vorticity. Specifically, we consider the 
cases where the vorticity is a 5-function (a point vortex) , or has small compact support (a 
vortex patch). Using a global bifurcation theoretic argument, we construct a continuum 
of finite-amplitude, finite-vorticity solutions for the periodic point vortex problem. For 
the non-periodic case, with either a vortex point or patch, we prove the existence of a 
continuum of small- amplitude, small- vorticity solutions. 

1. Introduction 

Consider the water wave problem of infinite depth modeled by the free boundary problem 
of the incompressible Euler equation 

(1.1a) 8fcv + (vV)v + Vp + 0e 2 = O, 

(1.1b) V-v = 

on a moving domain 

n t := {(x 1 ,x 2 ) £l 2 :i 2 <l + rj(t, x\)} 
for some profile function rj(t,xi), with the boundary conditions on the free surface dQt 

(1.2a) n t = -r] l vi+V2 

(1.2b) p = a 2 K, a>0 

where k = k(x\) is the mean curvature of the surface 

K ( x i) = r- 

(l + (r/) 2 ) 5 

Here — ge 2 = (0, —g) T is the constant gravitational acceleration. The kinematic boundary 
condition (1.2a) is the requirement that the normal velocity of the boundary agree with 
that of the fluid and condition (1.2b) means that we include the surface tension in our 
consideration. 



* The third author is funded in part by NSF DMS 1101423 and a part of the work was completed while 
visiting Courant Institute, New York University in 2011 and IMA, University of Minnesota in 2012. 
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Our goal in this paper is to seek traveling water waves in the form of 
(1.3) v = v(x\ — ct, x 2 ), rj = rj(xi - ct), = {(xi, x 2 ) G M 2 : x 2 < 1 + v( x i)} 
with compactly supported vorticity. We shall consider both cases 
(Loc) 7? and v decay at infinity and 

(Per) 7] and v are 2-7rL-periodic in x\, i.e. x\ S LS 1 and without loss of generality 



Since the arguments for both of these two cases follow a similar procedure, we will focus 
on the localized case (Loc) and then give outlines for the case x\ € S . 

The overwhelming majority of the research on water waves has been done in the irrota- 
tional setting, both for mathematical convenience and on physical grounds. The source of 
the convenience is clear enough: for velocity fields with gradient structure, the dynamics in 
the bulk of the fluid are captured simply by Laplace's equation. This allows the problem 
to be pushed to the boundary, where it can be recast in a number of ways as a nonlocal 
equation (e.g, as an integral equation following Nekrasov [15, 17], or using Dirichlet-to- 
Neumann operators as in Babenko [3].) Moreover, because Kelvin's circulation theorem 
states that an initially irrotational flow will remain irrotational absent external forcing, 
these considerations are physical. Today, the existence theory for steady water waves is 
very well-developed (cf., e.g., surveys in [12, 18, 21, 23]), though some important open 
problems still remain. 

Yet rotational waves are completely commonplace in nature. Indeed, the presence of 
wind forcing, temperature gradients, or even a slight heterogeneity in the density, generates 
vorticity. The study of rotational steady waves essentially begins with Dubreil-Jacotin in 
1934 (cf. [9]), but the entire topic lay dormant until the relatively recently. The main 
breakthrough came in 2004, when Constantin and Strauss developed a systematic existence 
theory for two-dimensional, periodic, finite-depth, traveling gravity waves (cf. [5]). Both 
the work of Constantin-Strauss and Dubreil-Jacotin begin with a simple observation: If 
there are no stagnation points in the flow, then one can use the stream function as a 
vertical coordinate to fix the domain. Doing so, one ultimately arrives at a quasilinear 
elliptic PDE on a strip (the quasilinearity coming from the change of coordinates), with 
a nonlinear boundary condition. Shear flows, i.e., those where the free surface and all 
of the streamlines are flat, can be easily described in the new coordinates, and so one 
can build small-amplitude solutions by a perturbative argument. Then, using a degree 
theoretic continuation method, Constantin and Strauss were able to obtain finite amplitude 
solutions. Following the publication of [5], many authors have been able to generalize this 
approach, and we now have a bounty of analogous results for other physical regimes (cf., 
e.g., [25, 13, 29, 30]), for a class of weak solutions (cf. [6]), and even for some types of 
waves with stagnation points (cf. [26, 10, 11]); a recent account of the field is given in [21]. 

Though they differ in many essential details, all of these works rely on the fact that 
the vorticity is constant along streamlines for steady flows away from stagnation. Indeed, 
even the papers dealing with the stagnant case, where this does not follow from physical 
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principles, impose it a priori. Another common feature is that they all use shear flows as 
their class of trivial solutions for the bifurcation argument. Taken together, these two facets 
make it impossible to construct waves with localized vorticity. To see why conceptually, 
simply note that in a shear flow, every point sits on a unique streamline extending to 
±00 in the direction of wave propagation. Along that streamline, the vorticity is constant, 
and hence it cannot be localized unless it happens to vanish identically. Naturally, small 
amplitude perturbations of these flows will share this property. Indeed, if there are no 
stagnation points, all of the streamlines are unbounded. With stagnation, it is possible 
that some of the streamlines become closed upon bifurcation — a phenomenon referred to 
as cats' eyes — but some unbounded streamlines will necessarily persist. 

Thus a vast gulf exists between the well-studied irrotational steady waves, where of 
course the vorticity vanishes identically, and the current literature on rotational waves, 
where it cannot be allowed to vanish at infinity. The objective of the present paper is 
to address this gap, and, in a sense, our approach weds the two outlined above. We will 
be able to describe traveling waves where point vortices or eddies are suspended in an 
otherwise irrotational flow. In particular, the traveling waves with compactly supported 
vortex patch found in this paper have finite energy. 

The vorticity u of a 2-d velocity field v is defined to be the distribution 

(1.5) u> := d Xl v 2 - d X2 vi. 

When a; is a finite measure, we may define the vortex strength to be e = u>(£l), and in 
particular, 



In the interior of the fluid, the incompressible Euler equation can be expressed in its 
vorticity formulation which takes the following form for traveling waves 



It means that the vorticity is transported by the flow. We consider two problems 

(PtV) oj{x) = e5(x), a point mass away from the fluid boundary, without loss of generality, 

whose concentration point is taken as the origin. 
(VPa) D : = supp (oj) CC and oj is smooth on D which is near the origin. 

In fact, we consider both the localized (Loc) and periodic (Per) cases for the point vortex 
problem (PtV) and only the localized case (LoC) for the vortex patch problem (VPa). The 
periodic case (Per) can also be considered for the vortex patch problem (VPa), but the 
computation is too tedious and we simply skip it in this manuscript. Vorticity equation 
(1.7) can be interpreted in the distribution sense when oj £ L\ oc as in the case of vortex 
patch (VPa). 

If for a solution of the incompressible Euler equation oj(t, x) = e8^u\(x) is a point mass 
concentrated at x(t) for some t, the general principle that the vorticity is only transported 
following the fluid flow suggests that oj might remain as such a point vortex for all t. The 
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question is, as the singular vorticity generates a singularly rotational part 

T7n(X2 - X 2 ,Xi - Xi) 

\x — x\ z 

of v, following what vector field should the vortex point move? Since the above singular 
vector field is purely rotational and does not move that particle at x(t) away, it is reasonable 
to expect that the dynamics of the vortex point is governed only by the remaining smooth 
part of v i.e. 

d t x = (v - . J _. 2 (x 2 - x 2 ,x 1 - xi) T )\ x=x . 

This well known result is rigorously established by considering a family of solutions whose 
initial vorticity limiting (weakly) to a 5-function (cf. [16, Theorem 4.1, 4.2]). For traveling 
waves where x(t) = cte%, this translates to the following weaker form of (1.7) 

(1-8) cei = (v - j^(-x 2 ,xi) T )\ x=0 - 

Indeed, our own analysis of the vortex patch confirms this for we found that the traveling 
speed c for the travel vortex patch satisfies (1.8) as the diameter of the patch converges to 
zero. 

As another indication how traveling waves with point vortex or vortex patch are different 
from those close to shear flows, it is worth pointing out that (1.8) and (1.7) imply that the 
vortex point or some point in the vortex patch is a stagnation point. That is, the physical 
horizontal speed coincides with the traveling speed c. Usually stagnation only occurs as a 
limiting case in traveling waves constructed via bifurcation from shear flows. They often 
coincide with the development of a singularity. 

Note that the water wave problem is invariant with respect to reflection in x±. For 
simplicity in this paper, we consider symmetric traveling wave solutions (fi, v) where f2 
and v\ are even in x\ and v 2 is odd in x\. 

Our main theorems are outlined in the statement below, while their more precise versions 
are given in Section 2. 

Main Theorem. (1) For < |e| <C 1, there exists a unique traveling water wave 
solution which is even in x\ with small amplitude and small velocity whose vorticity 
is given by a delta mass of strength e away from the surface. 

(2) This solution curve of traveling water waves with a point vortex can be extended 
globally with one of the possibilities: 

• either ij or v becomes unbounded, 

• the point vortex location becomes arbitrarily close to the water surface along 
the solution curve, or 

• there exists a nontrivial irrotational traveling water wave, with gravity and 
surface tension, with an interior stagnation point. 

(3) For < |e| <C 1, there exist traveling water waves which is even in x\ with small 
amplitude and small velocity whose vorticity of total strength e is compactly sup- 
ported in a small disk-like region away from the water's surface. 
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The construction of small traveling water waves with a point vortex, which is given in 
section 3, is based on a fairly simple implicit function theory argument. 

In section 5, we use a degree theoretic global bifurcation argument to extend the local 
bifurcation curve beyond the neighborhood of and into the finite-amplitude, speed, and 
vorticity regime. Part (2) is the result of this process, an alternative theorem in the spirit 
of Rabinowitz [19]. The global bifurcation curve must either be unbounded, with the 
separation between the point vortex and the surface limiting to along some sequence, 
or contain a nontrivial irrotational traveling water wave with gravity and surface tension, 
whose velocity at some point is equal to the wave speed c. Traditionally with global 
bifurcation arguments, one begins with a statement of this type, and then uses a priori 
estimates and nodal arguments to rule out one or more of the alternatives. This was the 
procedure used by Amick, Fraenkel, and Toland to prove the famous Stokes' conjecture (cf. 
[2]), and it was the means by which Constantin and Strauss showed that the limiting waves 
of their global continuum must approach stagnation (cf. [ ]). But, both of these results 
are for waves without surface tension, which is important because it renders the maximum 
principle dependent nodal arguments tractable. For rotational capillary waves, one of the 
authors has proved some theorems that go beyond simply the Rabinowitz- type (cf. [28]), 
but the large-amplitude regime is still mostly open. Intuitively, though, point vortices are 
nearly irrotational so that one might hope that Theorem 2.2 can be extended to match the 
state-of-the-art for irrotational waves. This is a very interesting question, but beyond the 
scope of the current paper. 

The third part of the main theorem, proved in Section 4, deals with the other class of 
vorticity that we study, the vortex patch. These are solutions where uj £ L 2 (0) n L l (Cl) 
and is supported in a compact region D CC fi. One can view the point vortex as a limit 
of vortex patches when the size of the patch taken to zero. We will require that uj be 
continuous on f2, smooth in D, but not necessarily across dD. Reformulating the vorticity 
equation in terms of the relative stream function, what results is a nonlinear elliptic free 
boundary problem, very much in the same vein as Constantin-Strauss. However, rather 
than perturbing from a shear flow, our point of bifurcation will be radial solutions on a 
ball. Moreover, we use conformal mappings to fix the support of the patch rather than 
streamline coordinates. By construction, then, the streamlines on which the vorticity is 
nonzero will be closed. On the other hand, the boundary motion of the air-water interface 
is dictated by (2.5) and (2.8), just as in the irrotational setting. To couple the two requires 
that some matching be done on dD. 

A few remarks about the hypotheses and possible extensions. For the vortex patch 
problem, the interplay between the interior dynamics and the boundary motion are more 
intricate, which is why we need the optimal regularity furnished by fractional order Sobolev 
spaces. Concerning the evenness assumption, while it turns out to simplify the computa- 
tion, it should be noted that this type of symmetry is often expected for traveling waves. 
Indeed, it is known that in a number of regimes, all traveling waves (with monotonic 
profiles) have an axis of even symmetry (cf., e.g., [8, 24, 4, 24, 27]). 

A global continuation along the lines Theorem 2.2 seems an order of magnitude more 
difficult to carry out for the surface 5^\ oc . The main obstruction is that, as can be seen 
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in section 5, the global theory requires being able to prove compactness properties of the 
linearized operator not only at 0, but anywhere along the continuum. While this is not a 
problem for the free surface equations, the elliptic PDE for the vortex dynamics requires 
significantly more finesse. We therefore do not consider this issue in the present work. 

The structure of the paper is as follows. In Section 2, we start to reformulate the 
problems into forms more suitable for our analysis. In Section 3, we give the proof of local 
bifurcation for traveling waves with a point vortex. The global theory of the point vortex 
case is developed using a degree theoretic argument in Section 5. We address the vortex 
patch problem in Section 4. Additionally, some auxiliary lemmas and technical facts are 
collected in an appendix. Finally, though we shall always define a notation when it is 
introduced, we include a table outlining our conventions at the end of the appendix. 

2. Framework and main theorems 

In this paper, we seek traveling waves near the trivial state, namely |v| <C 1 and \rj\ <C 1 
where rj is given in the definition of f2 in (1.3). In the interior of f2, the incompressible 
Euler Equation of traveling waves is equivalent to the vorticity equation (1.7) (possibly in 
a weak sense). In addition, there are two more conditions on d£l, the Bernoulli condition, 
and the kinematic condition, which come from (1.2a) and (1.2b), respectively. We will 
introduce the stream function and the velocity potential to reformaute the problem. 

Stream functions. On the one hand, in two space dimensions, it is convenient to write 
a divergence free vector field v as the skew gradient of a stream function 

v = V ± ^=: (-d X2 4>,d Xl iP) T , A* = u. 

In the periodic case (Per), even though the domain Q, is homotopic to S , and thus not 
simply connected, the decay of v at x 2 = — oo and the periodicity of v\ ensure that \t 
is single valued. Based on the convolution with the Newtonian potential, the vorticity oj 
naturally generates a corresponding part eG(x) of the stream function where G satisfies 

AG = ^_«y(._2e 2 ) 

along with 2-7rL-periodicity in x\ in the periodic case (Per). More explicitly, in the localized 
case (Loc), 

(2.1) G = (^-logH)*u;-^log|--2e 2 |. 

In the periodic case (Per) we need to sum up the above stream functions generated by the 
vorticity in each period and thus for x\ £ [0, 2irL), 

(2.2) G = — P.V. V ((log | • -2A:7rLei|) * — — log | ■ -(2kirL, 2) T |) . 

The 'P.V.' above means the principle value which ensures certain convergence of the summa- 
tion. In the point vortex case (PtV) when u(x) = eS(x), obviously the above convolutions 
simply yield a logarithm. While note (2k7rL, 2) T ^ Q as rj <C 1 is assumed, the purpose of 
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the second logarithm term in the summand above is to ensure that VG has better decay 
at x = oo, namely VG = O(yjp) m an cases if uj is compactly supported. Physically, it 
corresponds to the standard trick of introducing a phantom point vortex in the air region 
with equal but opposite vortex strength to correct for the lack of integrability of the New- 
tonian potential in M 2 . With G given as above, there exists a harmonic function on O 
such that 

(2.3) v = V^f = V X (ipn + eG). 
We will seek traveling wave solutions with £ H 1 ^). 

Velocity potentials. On the other hand, an irrotational vector field can be (maybe 
locally) written as a gradient field. In the localized case (Loc) where f2 is simply connected, 
obviously any € H 1 ^) has a conjugate harmonic function € H 1 ^) such that 
V ± '*Ph = V^h. In the periodic case (Per) where x\ G S 1 , the decay assumption on Vif)n 
as X2 — > —go and the harmonicity of imply that the circulation ffV^ipH ' (is vanishes 
along any closed curve. Therefore, ipu also has a single valued conjugate harmonic function 

<ph e h 1 ^). 

Since the vorticity u) of the traveling wave solutions (v, if) under our consideration is 
supported away from the free surface 30, V^G is also irrotational in a neighborhood of 
c?0. Therefore, there exists a function O (multi-valued in the periodic case (Per)) such 
that VB = V-'-G and we have the decomposition 

(2.4) v = V$ =: Vipn + eVe. 

Bernoulli equation in terms of steam functions. The incompressible Euler equation 
of irrotational velocity fields and boundary condition (1.2b) imply the Bernoulli equation 
of traveling waves on d£l 

-cd xi <& + 7^|V$| 2 + gx 2 + c?k = b on {x 2 = 1 + r)(xi)} 

where b is a constant. Though the velocity potential $ might be multi- valued as pointed 
out above, note here only is present which is single- valued. This fact allows us to write 
the Bernoulli equation in terms of the stream function only 

c(d X2 ^ H + ed X2 G) + -\Vi/)n + eVG| 2 + gx 2 + a 2 K = b on {x 2 = 1 + v( x i)}- 
Let ift be the trace of V% on the free surface i.e., 

1p(xi) = Tpu(xi, 1 + 

and thus ?/% := T-L(r])ip is the harmonic extension of ip to the fluid domain (cf. Lemma 
A.l). From the Bernoulli condition, we arrive at the following equation of only the variable 
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X\\ 



(2.5) 



= c (g( V )^ + e(-r/, if ■ VG) + \ (g(r,)^ + e(-rf, 1) T ■ VG) 2 



1 



2(1 + (V) 2 ) 

where dG is evaluated at {x\, 1 + rj(xi)) and 



(V' - f/Sfaty + e(l + (V) 2 )^iG) z + g(v + 1) + a 2 /c(r/) - b, 



and M(rj) is the Dirichlet-to-Neumann operator on f2; we recapitulate some of the properties 
of these operators in Lemma A.l. The constant b can be determined explicitly. In the 
localized case (Loc), by taking x\ = oo, we obtain 

(2.6) b = g. 

In the periodic case (Per), (1.4) implies 

1 



b = b(e;rj,ip,c) =g + 



2vrL 



wL 
-L 



c (G(rj)ip + e(-V, if • VG) + i (S( ?? )V + e(-r/, 1) T • VG) 2 



1 



ty' - ^(77)^ + e(l + (t/) 2 )SxxG) 2 + a 2 



2(1 + (r/) 2 ) 

In the above integral, the k term integrates to due to the periodicity of rj. Moreover, 

/tyL r r 

[G(rf)il> + e(-rf, 1) T • VG] dxi = / N • V(^ + eG) ds = / udx = l. 



Therefore we obtain in the periodic case (Per) 



b = b(e;rj,ip,c) =g + ce + 



1 



(2.7) 



2vrL 



(e(# + e(-r,',l) r -VG)' 



2(1 + (m 



(V , -^(r ? )V + e(l + (r? / ) 2 )^ 1 G) 



dx\. 



Kinematic equation. In terms of ip and G, the kinematic condition (1.2a) for traveling 
waves is written as 



{21 



= crf + + e(l,r/) T ■ VG. 



Equations (2.5) and (2.8) in terms of the velocity are derived, e.g., in [22] for the un- 
steady problem, here we have simply adapted them to the steady regime and transformed 
it using the stream function which is more convenient when the interior vorticity is treated. 



Point vortex problem. In the localized case (Loc), oo = e5(x) and thus 
(2.9) G(x) 



1 1 , 

7T lo S \ x \ - 7T lo S \ x ~ 2e 2l- 

Z7T Z7T 
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As discussed in Section 1, (1.7) is replaced by (1.8). In the localized case (Loc) where 
x\ G R, according to (2.1) and the symmetry of v in xi, we obtain 

(2.10) c = -(d x M(Q)-^- 

In the periodic case (Per) where x\ G LS 1 , (2.2) implies instead 

(2.11) G(x) = — P.V. ( lo S \ x ~ 2kirLei\ - log \x - (2kirL, 2) T |) 

fcez 

and 

(2.12) c= _ (axA)(0) _A ^ __. 

fc = — OO 

In the point vortex case (PtV), we will seek solutions (77, ^, c) with 77 even in x\ and y odd 
in x x , so that (2.5), (2.8), (2.6), and (2.10) in the localized case (Loc), or (2.5), (2.8), (2.7), 
and (2.12) in the periodic case (Per), are satisfied. 

To state the main theorems in the point vortex case (PtV), we introduce the following 
spaces of even functions. For the localized case (Loc), let 

H*(R) := {/ G H k (R) : / is even in x x }, 

ijg (R) be the corresponding homogeneous space, and 

(2.13) X := fl*(R) x (H^(R) n H*(R)) x R. 
For the periodic problem (Per), we simply replace X with X per , 

(2.14) X per := H^LS 1 ) x fl* (LS 1 ) x R 
where 

H^LS 1 ) := {/ G tf^LS 1 ) : / has mean 0, even in 

Our first theorem establishes the existence and uniqueness of a curve bifurcating from the 
trivial solutions. 

Theorem 2.1 (Point vortex local bifurcation). Consider the traveling water wave problem 
with a point vortex at the origin (2.5), (2.8), (2.6), and (2.10) in the localized case (Loc), 
or (2.5), (2.8), (2.7), and (2.12) in the periodic case (Per). The following statements hold. 

(i) There exists eo > and a C°° -curve of solutions 

< gioc = {(e > f7(e),V(e),c(e)) : |e| < e } C R x X (orRxI pcr ) 
for any k > | with 

(0,r/(0),^(0),c(0)) =(0,0,0,0). 

Moreover, in a sufficiently small neighborhood of in R x X (or R x X per ), <^i oc 
comprises all solutions. 
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(ii) In the localized case (Loc), the solutions have the asymptotic form 



anc 



4vr 2Vy XlJ y {l+x 2 ) 2 ' 

(iii) In the periodic case (Per), the solutions have the asymptotic form 

c(e)=ec + o(e 2 ), \^(e)\ Hk = 0(e 3 ), \ V - e 2 V *\ Hk = 0(e 



where 



co 



1 00 1 
4vr ^ k 2 ir 2 L 2 + l' 



-oc 



r,*=-(g- cfdlX 1 (co^ 2 G + \{d X2 G) 2 - ^ J* cod X2 G + ±(d X2 G) 2 d Xl ), 

G is defined as in (2.11), and VG is evaluated at x<i = 1. 

Remark 2.1. (a) Note the above existence and uniqueness results hold for all k > |, so the 
solution (ji, ip) are actually C°° functions. 

(b) The above local uniqueness is stated in the framework of (77, -0, c) € X (or X per ). 
However, if a velocity field v satisfies V • v = 0, V x v = e5, and v £ L 2 outside any 
neighborhood of 0, it can be written in the form given in (2.3). Therefore, the local 
uniqueness holds in the space of such vector fields. 

The proof, which is given in section 3, is based on a fairly simple implicit function 
theory argument. In section 5, for the periodic case (Per), we use a degree theoretic global 
bifurcation argument to extend the curve beyond the neighborhood of and into the finite- 
amplitude, speed, and vorticity regime. The result is the following alternative theorem in 
the spirit of Rabinowitz [19]. 

Theorem 2.2 (Point vortex global bifurcation). For k > 3 in the definition of X pcr , there 
exists a connected set 'rf C R x X pcT of solutions to (2.5), (2.8), (2.7), and (2.12) with 
'rfloc C ' . One of the following alternatives must hold: 

(i) There is a sequence {(e n ,7] n ,ip n ,c n )} C "«f that is unbounded mix X per ; 

(ii) there exists a nontrivial irrotational (i.e. e = 0) traveling wave solution with stag- 
nation point at the original, or 

(iii) along some sequence in {(e n ,rj n ,ip n ,c n )} C , we have %(0) — > —1. 

Localized waves with vortex patch case. In this case, we will look for solutions 
(t],i/j,uj,c) of (2.5), (2.8), (2.6), and (1.7). In fact, (1.7) is essentially only required on 
D := supp (cj) while it should be satisfied on Q in the distributional sense. Since D CC 
is an unknown, through a more elaborated procedure using conformal mappings in Section 
4, (1.7) will be further transformed to a form more suitable for our analysis. Using a (local) 
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bifurcation theory argument, we are able to construct a curve of small-amplitude, small 
vorticity, slow speed, and small patch solutions in the neighborhood of the trivial solution. 

Theorem 2.3 (Vortex patch local bifurcation). For any s > 3/2 and integer ko > 1 
or ko = oo, there exists a three-dimensional C k ° surface of solutions to the vortex patch 
problem 

^loc = {('q(S,e,T),ip(S,e,r),u(S,e,T),c(d,e,T)) : (e,8,r) G U}, 
where U := [0, eo) x (0, 6q) x [0, to) for some eo, So, tq > 0, 

^loc C H S C +1 (R) x (H s e (R) n i?j(R)) x H^(R 2 ) x R. 
77ie parameterization is such that 

(77(0,0,0), ^(0,0,0), w(0,0,0),c(0,0,0)) = (0,0,0,0), c = e (__L + o( e + ( y)). 

Moreover, oj G H S+ ?{D) where D := supp (cj) C and is an £P closed curve with 
asymptotic form: 

dD(e,5,r) = {5 (cos6 + Tsm(26),sm6 -t cos (26)) + 0(5 2 {5 + e)) : G [0,2vr]}. 

Remark 2.2. 1.) Here we see that, for fixed e > 0, as 5 — > 0+, the wave speed c converges 
to the wave speed of that of the point vortex problem. 

2.) In fact, such a smooth family of solution is found for any fixed vorticity strength 
function 7 from a large class of functions to be introduced in Section 4. By choosing 
different 7, potentially a large family of such traveling waves can be found which differ in 
the 0(<5 2 (<5 + e)) term in the above. 

3. Small amplitude waves with a point vortex 

Let us first turn our attention to the point vortex problem. We will only focus on the 
localized case (Loc) and the proof for the periodic case (Per) follows from exactly the same 
procedure. Since the traveling wave solutions we are seeking are bifurcated from the trivial 
solution, r\ and tp are of order 0(e). In principle, the wave speed c can be anything for the 
trivial solution, equation (2.6) implies that c is also of this order, which is one of the major 
differences between the problems with localized vorticity and the one near the shear flows. 
Therefore it is more convenient to rescale 

rj = efj, ip = eip, c = ec. 

Using abstract formalism, we can express the governing equations (2.5), (2.8), (2.6), and 
(2.10) as 

F(e; fj, -ip, c) = 0, 



12 SHATAH, WALSH, AND ZENG 

where T = (Ji, T2, Tz) ■ K x X — > Y, and 

T x := ce (g(efj)i> + if -Vg) + ^ (g(e^ + (-er?', if ■ Vg)' 

? - efj'g(efj)tp + (1 + (er?') 2 )a xi G) 2 + 577 + —K(efj), 



(3.1) 2(1 + (ef/') 2 ) 

F 2 :=cefj' + + {l,efj'f ■ VG 

jP 3 :=c+(^ 2 ^)(0) + i- 

where G is defined in (2.9), VG is evaluated at x 2 = 1 + efj, b is replaced by g as in (2.6), 
the space X is defined as in (2.13); and Y is taken to be 

Y := H*- 2 (R) x (iJ^-^M) n # e ~^(M)) x R. 

Obviously small traveling waves with a point vortex at the origin corresponds to a zero 
point of J-. We are now ready to prove the local bifurcation theorem for point vortices. 

Proof of Theorem 2.1. Recall the VG term has the explicit form 

Thus the VG term depends analytically on {fj,ip,c) for |e| <C 1. Similarly, the operators 
T~L, G, and k are also smooth (cf. Lemma A.l), implying that J- is of C°° class (in fact they 
are analytic, but this is beyond our needs.) 

For e = 0, is clear that the following point is in the zero-set of T: 

770 = 0, V'o = -G| X2= i = 0, co = -— . 

Now, denoting 

{DP) := (Dfj, , D 5 )T(0; 0, 0, - -L ) , 

a simple computation reveals that 

/ g - a 2 d^ \ 

(DF) Q = a ai G£(X;F). 

V (d X2 (H(0),-))\ m 1) 

Each diagonal entry here is invertible (see Lemma A.l), and thus Z).F(0) is an isomorphism. 
The implicit function theorem immediately implies the local existence and uniqueness of a 
curve of zero points of T parametrized by |e| <C 1. Direct expansion shows c = — 4- + 0(e), 
$ = 0( e 2 ) and 

r/ =6( 5 - a 2 ^)"^ - co^G - \{d X2 G? + \{d xl Gf) + 0(e 2 ) 
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where VG was evaluated at X2 = 1. This in turn yields c = — 4- + 0(e 2 ). The proof for 
the localized case (Loc) is complete. □ 

4. Small amplitude waves with a vortex patch 

We now consider traveling waves in the localized case (Loc) where x\ £ R and the 
vorticity ui is supported in a compact region D := supp (w) CC 0. As outlined in section 
2, this problem is equivalent to the system given by (2.5), (2.8), (2.6), and (1.7) with 
unknowns (rj,tp,u,c) with e := J n ojdx as a small parameter. Since ui is supported on an 
unknown domain D, we first further rewrite the problem to address this difficulty 

4.1. Reformulation. As one of the unknowns, we consider the vortical region D as being 
a perturbation of B$, the ball of radius 5 centered at the origin, where 5 is fixed and small. 
As we assumed d£l = {x2 = Tj(xi) + 1} and \n\ « 1, there is a positive separation between 
the patch D and the air-water interface. 

As given in (2.3), recall VP = ipu + G is the stream function of v where G is given in 
(2.1). Define the relative stream function 

/ := ^ + ct 2 , 

then the vorticity equation (1.7) takes the equivalent form 
(4.1) V/-VA/ = V ± /-Vw = 0. 

Before proceeding further, it will be convenient to scale the physical variables, the relative 
stream function, and the vorticity. With that in mind, define Dq := and /, uj : Dq — > R 
by 1 

(4-2) f(x) =: 6/(|), u(x) = J^(|). 

Note L>o is thus a small perturbation of B\ and that, as Af = uj, 



Af = uj in Dq, / Af dx = I ujdx = l. 

J D J D 

One way to ensure that the vorticity equation (4.1) is satisfied is to require that the vortex 
lines and streamline coincide, i.e. 

A/ = £ = 7 (/) in D , 

for some function 7 called the vorticity strength function. As Dq = supp(S), 8Dq is a 
vortex line (lj\q£i = 0) and thus a streamline where / is a constant. Without loss we may 
assume that f\dD = so that, in total, the scaled relative stream function / € Hq(Dq) is 
the solution to the elliptic PDE 



(4.3) 




in D 
on 8Dq. 



^Here and elsewhere in this section, we adopt the convention that quantities scaled in 8 and e are 
identified with tildes. 
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Since the above elliptic boundary value problem can be approximated by the one with 
Dq = Bi, we impose some conditions on 7 so that it has a non-degenerate solution when 
Dq = B\. Namely, we assume 

(4.4a) 7GC JV (R), 7(0) = 0, 7 ; (0) < 0, 7 > on R", 

(4.4b) there exists a negative radial solution /* to (4.3) with Dq = Bi, 

and 

(4.4c) A — 7'(/*) is non-degenerate, 

where the integer iV > s + ko + 9 with s and ko being given in Theorem 2.3. We will look 
for traveling wave solutions with Dq close to B\ for each such fixed 7. 

Solving (4.3) will allow us to determine the vorticity, but before we can do that, we must 
represent the domain Dq in an analytic way, which is accomplished by using near identity 
conformal mappings between B\ and Dq. Let T be a conformal mapping with domain the 
unit ball B± C C, and satisfying 

(4.5a) d 2 T = in B u T(0) = 0, r'(0) = 1. 

By identifying z = x\ + 1x2 £ C with the point (xi,X2) £ R 2 , we may view the dilated 
domain Dq as the image of B\ under T, Dq := T(Bi), and the unsealed domain D := 5Dq. 
As we are interested in vortex patches that are perturbations of B$, we think of \T(z) — z\ 
as being small throughout Do. 

Let us now briefly motivate (4.5a). The first statement is just the conformality, while 
the second fixes the origin. The third is made in order to eliminate a certain redundancy. 
Observe that for each a > 0, (S,T) and (5/a,aT) each result in the same patch D. By 
fixing r'(0) = 1, we exclude all but a = 1. 

In addition, since we look for traveling waves with the fluid domain and the stream 
functions even in x\, we require that D is symmetric over the xi-axis, which, stated in 
terms of T, is equivalent to 

(4.5b) ReT is odd in x±, and ImT is even in x\. 

It is clear that there is a one-to-one correspondence between a symmetric domain Dq close 
to B\ and such a near identity conformal mapping V satisfying (4.5a) and (4.5b). 

The connection between (4.5b) and the symmetry of the domain can be seen as follows. 
As T is analytic, we may express it as a power series 

00 

(4.6) T(z) = Y,a n z n = z + ^a n z n . 

n=l n>2 

That ao = and a\ = 1 follows from (4.5a). Consider the term a n z n , for some n > 1. 
Denote 6 n := arga n and 9 := argz, which implies 

arg(a n z n ) = nd + n . 
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Symmetry of D over the xi-axis translates to the requirement that 

cos (n(7r — 6) + n ) = — cos (nO + 9 n ), 
sin (n(7r — 6) + n ) = sin (n# + n ). 

Expanding the first of these identities yields 

(-l) n+1 cos {n6 - 6 n ) = cos (nO + 9 n ). 

From this it is apparent that cos 6 n = for n even, cos 6 n = ±1, for n odd. That is 

(4.7) Rea 2 „_i = a 2n -i, Rea 2n = 0, for n > 1. 

As a consequence, we have (4.5b). Moreover, the above property of a n implies that T takes 
the form 

oo 

(4.8) T(z) = z + Y, * n -%^\ Pn e R. 

n=2 

Note that the correspondence between a conformal mapping T defined on B\ and the 
real part of its trace on 8 1 is one-to-one and onto. Let 



oo 

(4.9) P = (3(e)=Re[T(e ie )-e ie } = ^2p n cos{-(n-l)Tr + ne), 8 £ S 1 . 

n=2 

the trivial solution corresponds to (3 = 0. For each s > 0, define the space 

oo 

(4.10) X s := {/? G ff'CS 1 ) : /? = £ /3„ cos (-(n - 1)tt + nfl), {/3 n }~ 2 C M}. 

n=2 

Obviously /3 G ff s if and only if ^~ =2 ^ 2s /?n < oo. 

We consider f3 G H S (S 1 ) to be the unknown describing the shape of the vortex patch, 
and write V = T((3) to emphasize that V is defined in terms of (3. It is not hard to show 
that r(/3) is smooth and bounded in the appropriate spaces, cf. Lemma B.l. Where there 
is no risk of confusion, we shall abuse notation and suppress this dependence entirely. 

Letting F := f o r(/3), the semi-linear problem for / becomes 

AF = |^r| 2 7 (F) on B x 
F = onS 1 . 

We have thus managed to fix the domain. Note, however, that the above equation will not 
suffice as a definition of F, because we need in addition that f B ^ AF dx = J Do Afdx = 1. 

With that in mind, we instead consider the problem of (F, a) 



(4.11) 



AF = a|^r| 2 7 (iF) on B x 
a 

F = 0^ on S 1 

f Bi AFdx = 1, a > 0. 
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Now we proceed to transform the system given by (2.5), (2.8), (2.6), and (1.7) with 
unknowns (rj, tp, u, c) to a system of equations with unknowns (rj, ip, /3, c). The key here is 
(1.7) which is already expressed by the above semilinear elliptic problem. 

Given /? S X s , by the nondegenercy assumptions on 7 in (4.4), for f3 G X s in a neigh- 
borhood of 0, there exists a locally unique solution (F,a) to (4.11). Indeed, we show in 
Lemma B.l that the solution will depend smoothly on f3. Once T and F are known, the 
scaled vorticity support Dq = T(B\) and define u through 

(4.12) uoT = \d z T\~ 2 AF in D and u = in (D ) c 

which is expected to be the scaled vorticity. Accordingly, we define the mapping 

(4.13) h : p e x s ^ [ log |r(-) - t(zM-f(z))\t' (z)\ 2 dz 

2vr J Bl a 

where F is the unique solution to (4.11). The smoothness of H is established in Lemma 
B.3. 

Chasing definitions, H(/3) is expected to be the scaled trace on ODq of the rotational 
part of the relative stream function corresponding to the vortex patch. We will need to 
choose r appropriately so that we can take the solution F to the semilinear problem and 
find a corresponding relative stream function defined globally in f2 for the traveling wave 
solution we are seeking, for this, understanding the dependence of H on (5 is critical. In 
fact, given (77, ip, f3, c), define 



(4.14) 



where 



~ 1 c 
f(x) = -ty(Sx) + -5x 2 

= ipu(5x) + — (log I • I * uj) (x) — — log \6x - 2e 2 | + 80x2 + /J, 

ifn := -tpn, c:=-c, 
e e 



and n is chosen so that the mean of / o T on S 1 is zero. / is expected to be the relative 
stream function. Based on these definitions, we have 



and thus 



A/ = £5, #(/3) = ^(log|-|*£)or 



f(x) = ifj n (5x) + {H{fi) o r- x )(x) - — log \Sx - 2e 2 | + 5cx 2 + fi. 

Z7T 



If / o r coincide with F, which solves the elliptic problem (4.11), then clearly / satisfies 
(4.3) and thus it yields a velocity field v satisfying (1.7). Since AF = A(/ o T), F\§i = 0, 
and J gl / o T ds = 0, we only need 



= d e 



$ n o {5T({3)) +HImr(/3) - i- log \5T((3) - 2e 2 | + H{/3) 
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Note that the bracketed function, when considered on S 1 , is even in x\ and thus the right 
side of the above equation is odd in x\. This means that it enjoys the same symmetry as 
P does, but may include the n = 1 mode. Define 



(4.15) 



oo 1 

Y s := {p G H^S 1 ) : /? = /3 n cos (-(n - 1)tt + n9), 



{PnYr 



C 



n=l 



Let us now put this in the framework of the implicit function theorem. As in the point 
point vortex case, we consider the scaled unknowns 



1 



1 



V 



:V, 



„ c 
c = - 
e e 

as well as j3 which gives T, Dq = r(£?i), F, uj. Based on the above discussion along with 
the Bernoulli equation (2.5) and the kinematic boundary condition (2.8), the problem of 
traveling water waves with a vortex patch is transformed to a system in the form of 



where 



(4.16) 



srf, if • VG) + | (g{en)i> + (-erf, if • VG 
j'-erj , g(e?))$+(l + (erj') 2 )d xl G 



T 2 



2(1 + (e^') 2 ) 
cerf + tfj' + (1, erj'f • VG 



2 „ a 2 
+ 9V H K(erj), 



cSluiT(P) - —log\5T(P) 



2e, 



H(P) 



Here VG is evaluated at x 2 = 1 + er](xi), G is defined as in (2.1) with uj given by (4.2), 
(4.12), and (4.11), which are solely determined by p. Any zero point of J- defines / by 
(4.14) which after rescaling (4.2) yields the velocity field v = V _L (/ — cx 2 ) € L 2 (£l), where 
£1 is given by erj, for the traveling wave solution. To find zero points of J 7 , we will prove in 
the appendix (cf. Lemma B.4.) that H is a C ko+1 mapping from X s to Y s and T is C ko+1 
from M. 2 x X to y with parameters e and 5 where 

X := H° +1 (R) x (H*(R) n # C 1/2 (M)) x X s x R, 

^ 4 ' 17 ^ y := iJ e s_1 (R) x {H^\R) n ^~ 1/2 (1R)) x Y s ~\ 



4.2. Proof of local bifurcation for the vortex patch. The remainder of the section is 
devoted to constructing elements of the zero-set of J- with < e, 5 <C 1; these correspond to 
solutions to the vortex patch problem with nontrivial total vorticity. As a preliminary step, 
let us note that the first two components of J- are formally very similar to the corresponding 
operator for the point vortex case. Thus ^1,2(0, 5; 0, ip, P, c) = and the first two rows of 
the operator matrix DJ-(0,8;0,iJj, P,c) have invertible entries on the diagonal, and zero 
entries elsewhere. Moreover, when f3 = 0, our assumptions on 7 implies that there exists 
a unique, radial, negative solution (F*,a*) to the unperturbed problem found by taking 
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r = tBn the identity function on B\, in (4.11). In particular, the radial symmetry of F* 
immediately implies H(0) = and thus J-(0, 0; 0, 0, 0, c) = for any c. Summarizing these 
properties we have 

Lemma 4.1. Let J 7 be defined as in (4.16). Then .7-1,2(0, 0; 0, o, 0,c) = and 

DJi (0,05 0,0,0,2) = (5-a 2 ^, 0, 0,0), RF 2 (0,0;0,0,0,B) = (0, 0(0), 0, 0), 
where D = (D?j, D^, Dp, Dc). Consequently, at (0, 0; 0, 0, 0, c), 

A^i : #| +1 (R) -> H"| _1 (M), 

and 

z^F 2 : ff c s («) n # C 1/2 (M) -> iir 1 (R) n #~ 1/2 (M) 

are invertible. 

Because D^T\ and DrJ-2 are isomorphisms at 0, the main obstacle is in showing that 

DpFz = Dp(dgH(/3)) „_ Q (at <5 = 0) is an isomorphism. Before proceeding with that 
analysis, though, we pause to make some useful observations. 

Recall that we are seeking solutions where the vortical region is a perturbation of the 
ball of radius 5 > centered at the origin. Because F* is radial, and the linearized problem 
around it is non-degenerate, we will see shortly that it is helpful to consider the inverse 
of the linearized operator in polar coordinates. Assumption (4.4), the maximum principle, 
and the radial symmetry imply 

(4.18) d r F*{0) = 0, d r F* > for r G (0, 1]. 
Moreover, by the definition of total vorticity and our normalization, 

(4.19) 1 = f u* dx = 2vr f rAF* dr = 2vr f d r (rd r F*) dr = 2-ird r F*(l). 

JBj Jo Jo 

We define for each n > 1, 

(4.20) D\ n (r) := 2a*(c\ 2 + U r - ^ - 1 \^)y\ 1 {±-F*) r n ]. 
Proposition 4.1. Let H be defined as in (4.13). Let f3 G X s be given by 

n>2 

Then for z = e ie G S 1 , 

(4.21) [D P H{0)$] ( z ) = 1 Lj2l 1 - tt— ^(lH/Wi cos (^ + n6), 

11 n>l 71 

where 9\ n are real-valued functions defined in (4.20). 
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Proof. Let {/3a}a>o be a family in X s smoothly parameterized by A such that /3q = 0; 
let -{Ta}a>o be the corresponding family of conformal mappings uniquely determined by 
(4.9). Notice that this implies T = iBi(o) • Take {(F\,a\)} to be a smoothly parameterized 
family of solutions to (4.11) for T\, where (Fo,ao) = (F*,a*) and put 

h x :=H(p x ), J\ := F\ o T x . 

Chasing definitions, this means that F\ satisfies 

1 



A 



-i 



a x \d z T x \^(—F x ) 
ax 



Here A 1 indicates the inverse Laplacian with homogeneous Dirichlet conditions on -Bi(O). 
Let 



P := D X \ X=0 (3 X , 
For each z G S , 

h{z) = [D> 



T := D 



A|A=0 1 A> 



F := D 



A|A=0 



Fx, 



h := D 



hx. 



1 log \T x (z) - T x (z')\a x \T' x (z')\h(-Fx(z')) dz'] 



2tt 



ax 



X\X=0 n X 



A=0" 



Evaluating the derivative on the right-hand side reveals 

h(z) = — [ \og\z- z'\AF(z')dz' 
2? r Jb x 



+ 



1 

2^ 



Re 



(z-z>){T{z)-T{z>)) 



B, 



\z — z 



f\2 



AF*(z')dz' 



(4.22) 



i- f ( log |* - z'l AF(z') + Re T } Z '\ u 



By 



Z — Z' 



1 

2^ 



(I + 11). 



Likewise, the variation F solves the linearized problem 



(4.23) 



AF = i{^F*)F + 2a* Re [t'M^F* 

+a[ 7 (±F*) - iy(l^)^] 
a* a* a* 



/ Bi AFcix = 
F = 



in Si 

in Bi 
on <9I?i. 

To analyze (4.22) we make use of the analyticity of T x to expand it as a power series, 

T x = z + Y,a n (X)z n . 

n>2 
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The coefficients depend smoothly on A, and thus 

n>2 

Inserting this ansatz into the equation satisfied by F (4.23), we find 

AF = i{\F*)F + 2a* Re na n z n ~ 1 ]j(-^F*) 
a* ^— ' a* 

(4.24) n ^ 2 

+ a[ 1 (\F*)-±~ / '(\F*)F*]. 
a* a* a* 

For simplicity, let us reindex the expansion by taking b n := a n +i, so that 

(4.25) t = J2 b n z n+ \ F = + l)b n z n . 

n>l n>l 

Comparing (4.25) to (4.8) and (4.9), it is clear 

(4.26) b n = i n $ n+1 . 

Consider the second term on the right-hand side of (4.22), II, and use the radial sym- 
metry of oj* = a*^/(F* /a*) and the fact that J B uj*dx = 1, we obtain 

11=/ Re ~ y ^ ]u*(z') dz' = ! Re[Y J bnY.z j (z') n ~ j ]u*(z')dz' 
Jb 1 z-z J Bl ^ . =Q 

= Re £M / aS"V)dz' = Re— = Y p n+lC os(— +nd). 
z Jb, z tTi 2 



(4.27) 



z Jb i „>1 
Next we treat I, the first term in the integrand on the right-hand side of (4.22). 

1= lim / log | ^ - z'\AF(z')dz' 

r 0^l- J Bro 

= lim (f log\z- z'\V N F(z')da(z') - [ V N {\og\z - z'\)F{z') do{z'\ 

r -H- \ JdB ro JdBr 

where we integrated in parts twice and Vat denotes the outward normal derivative. It is 
easy to show that the second integral vanishes in the limit tq — > 1. Since we are integrating 
over a ball, Vtv = d r , and hence the polar representation of the integral for z' = r^e 1 ® is 

(4.28) 1= lim I(r ) := lim (r / log \e ld - r e l6 ' \(d r F)(r e i6 ) d0' ) . 

ro-^l- ro-»l- V Jq J 

Another way to write (4.23) is 



F = (A- 1 '(±F*))- 1 
a* 



2a* Re [f'M^*) + *(^ F *) ~ \f*J{\f*)) 
a* a* a* a* 
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where the inverse exists by our assumptions on 7, and is taken with homogeneous Dirichlet 
boundary conditions. As many of the quantities involved are radial, it is advantageous to 
express everything in polar coordinates. Write 



F = d + G 2 



where 



G 1 := (A- 7 '(1f*))" 1 



2a* Re [t'] 7 (-F* 



and G2 = F — G\ is radial. Writing the inverse in polar coordinates and expanding T' and 
using (4.25) and (4.26), we find that 



G!(re») = 2a* £> r 2 + \d r - £ - j'(±F*)) 



1 



*\\-i 



n>l 



7 (^F)r n Re((n + l)6, 



JHn(r)(n + l)/3 n+ i cos + n0). 



n>l 



Thus, for z' £ <9-B ro with argz = 9', 



{drG^z') = Y, K(ro)(n + l)/3„+i cos + n 0>). 



n>l 



Inserting this expression into (4.28), we find that 



n>l 



2k 



Y K(ro)r (n + l)$ n+ i / log |1 - r e^'-^\ cos (^ + nC) d& 



,nir 



+ r d r G 2 (r ) / logll-roe^'-^l^' 
J 

^ <(r )r„(rc + l)/3 n+ i / * Re []T -r^ e '" e ) cos + ^' 

n>l 7>1 J 



^±±rt n (r )$ n+1 r% +1 cos (™ + nO). 

ri>\ 



E 



Taking ro — > 1, this shows that 



(4.29) 



I = -7T 



E 

n>l 



77-1-1 ■ 777T 

_^3^(l)/3 n+1 cos(— +n0) 
n 2 



Finally, combining (4.22), (4.27), (4.28), and (4.29), we arrive at (4.21). This completes 
the proof. □ 

Lemma 4.2 (Fredholm map). Let H be defined as in (4.13). Then the following statements 
are true. 

(a) jV{DpH{§)) is one- dimensional and generated by the map 6 1— > sin (29). 

(b) &(deDpH(0)) = (# H> cos^)^ is codimension 1 in Y 5 " 1 . 
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(c) dgDpH(0) G L(X S ,Y S ) is a Fredholm operator of index 0. 

Proof, (a) In Proposition 4.1 a formula was derived for the variations of H in (3. We now 
proceed to identify its null space, which will require a closer look at d\ n . Adopting the 
same notation as in the previous lemma, let us additionally define for n > 1 



d 'r + \ d T 



-1 



7'(^n ) h&*)r n ] = i*n(r). 



n 

r 2 v a* ' j L ' x a* ' 2 



Explicitly, q n solves 
1 

(4.30) 



1 ~ r? 2 1 ~ 

d r (rd r q n ) - W(-F*) + -)q n = a * 7 (-F*)r" on (0, 1) 
r a* r z a* 



<7n(0) = g„(l) = 0. 

While g n (l) = comes from the definition of the inverse, the vanishing of q n at the origin is 
required in order for q n to be smooth at r = 0. These conditions imply q n G C N . Moreover, 
a blow-up analysis at r = 0, an ODE regular singular point, reveals q' n (0) = for n > 2. 

Claim I: q n is non-positive on (0,1). In particular, there is no interval (n,r2) C (0,1) 
such that 



(4.31) 



q n > on (ri,r 2 ), q n (ri) = q n (r 2 ) = 0. 



Seeking a contradiction, suppose that there exists < r\ < i'2 < 1 as above. Then from 
(4.30) we compute 



'■_> 



(4.32) 



''2 



'1 



1 ~ n 2 
d r {rd r q n ) - r( 7 '(-F*) + -)q n ) 
a* r z 



d r F* dr 



1-2 



u 1 



r{d z r F*)d r q n + rq n (d r F*)(j( — F*) + 



+ (rd r F*d r q n 



ri 



u 1 ^ 



drird'F*) - ridrF^'i-F*) + 



+ (rd r q n d r F*) 



dr 



Here the fact that q n vanishes at the endpoints has been used in the second integration by 
parts to cancel the boundary terms. Recalling that F* is radial and satisfies (4.11) with 
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r = , we see that 

d r (rd 2 F*) = d r (ra*j(—F*) - d r F*) = a*j(—F*) + r 7 '(— F*)8 r F* - 3 2 F* 
a* a* a* 

= -d r (rd r F*) + n '(—F*)d r F* - 8 2 F* 
r a* 

= -d r F*+rV(—F*)d r F*. 
r a* 



Inserting this identity into (4.32), gives 

pra ~, . f r2 1-n 2 



j r n+1 u*d r F* dr = j q n d r F* dr + (rd r q n d r F* 



ri 



This is impossible: the integrand on the right-hand side is non-positive by (4.31) and 
(4.18), yet <9 r g n (r 2 ) < 0, and d r q n (ri) > so that the second term is non-positive. We 
have proved the first claim. 

Claim II: The map j3 : 8 \— > sin (28) is in the null space of H := DpH(Q). First observe 
that setting n = 0, r 2 = 1 and proceeding as in the proof of the first claim leads to the 
identity 

(4.33) C r n+1 £j*d r F* dr = C qJ—^-d r F* dr + (d r q n )(l)(d r F*)(l). 

Jo Jo r 

In particular, when n = 1, (4.33) implies 

d r F*(l)d r qi{l)= [ r 2 u*d r F*dr 
Jo 

r 1 i ~ 

= / r 2 {-d r (rd r F*))d r F* dr 
Jo r 
1 / ~ \ 2 

= 2<V*(ir 



Since d r F*(l) ^ 0, we may divide to find 

(4.34) d rqi (l) = \d r F*(l). 

Combining (4.34) and (4.19) gives simply that 

m> 1 (l)=2q' 1 (l)=d r F*(l) = ±-. 

Now, if (3(8) = sin (28), then the coefficients are given by /3 2 = — 1 and f3 n = 0, n ^ 2, and 
thus from (4.21) we compute 

H(8) = - — [1- 27rJKi(l)]sin0, 
2tt 

which vanishes identically. Hence 8 i— > sin (28) E jV(D^H({))), as claimed. 
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Claim III: Let $ G X s be given by 

71 

/3:6>^cos(-7r + (n + l)60. 

Then /3 is not in the null space of DpH({)) for any n > 1. To prove this, consider again the 
identity (4.33). By the first claim, q n is non-positive on (0, 1), and so 



d r q n {l)d r F*{l) < f l r n+l u*d r F* dr. 
Jo 



Now, using the PDE satisfied by F* 



i 



r n+l u*d r F* dr = I r n+1 \d;F* + -d r F* ) d r F* dr 



1 



< 1 J (2r n+l d r F*d 2 r F* + {n+ l)r n (d r F*f^ dr 



= \{d r F*{i)f. 

Thus, 

d r q n (l) < \d r F*{l) = i-. 

As an immediate consequence, we have that 

n + 1„, , s n + 1 „ , . n — 1 1 

4.35 1>1-tt 9^(1) = 1 - 2tt d r <?„ 1 > — — > n>2. 

n n in 4 

Using the variation formula (4.21), we compute that 

DpH(0)p] {e w ) = i-[l - 2^(1)] cos + n0). 

This function cannot vanish identicaly, since the bracketed quantity is nonzero by (4.35). 
It follows that no $ of this form can be in the null space. This completes the proof of the 
third claim. 

Part (a) is a direct consequence of Claims II and III. Next consider (b). From (4.21) we 
see that for any $ G X s , 

(4.36) [doD H(O)$\ = i-^n/3 n+1 [7r^X(l) " 1] cos (^ir + n#). 

^ n>l n 

Part (c) follows immediately and the proof of the lemma is compete. □ 

We are now in a position to prove the main theorem for the vortex patch problem, 
namely the existence of small-amplitude solutions, via Lyapunov-Schmidt reduction. 

Proof of Theorem 2.3. Let 

Xq = {(O,O,fsin(20),c) :r,?GlR}, X x = X^ C X 

y = span{(0, 0, cos 6)}, y 1 =% i c3' 
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and Tlx and Tly be the associated projections to Xq and 3^o> respectively. From the Implicit 
Function Theorem and Lemmas 4.1 and 4.2, there exists a neighborhood IA of the line 
{(e = 0, 5 = 0;t = 0, c)} C R 4 and a (7 fco+1 mapping (Aft + 1 given in Theorem 2.3) 

u = (u±, U2,us,c) : U — > X\, u(0, 0, 0, c) = (0, 0, 0, c) 

such that in a neighborhood of the line {(e = 0, 5 = 0; rj = 0, tfj = 0, f3 = 0, c)} C R 2 x A' 

J 7 = iff rj = u±(e, 5, f, c), ^ = it2(e, <5, f, c), /3 = fsin(20) + 1x3 (e, 5, f,c). 

What is left to achieve is to find a solution of 

(4.37) = h(e,S,T,c) := ^- [ cos (6) T 3 (e, 8; u(e, 5, 6r,c) + (0, 0, <5rsin(20), 0)) d9. 

w Jo 

Note we have built the scaling f = St and \Fz into the definition of h. In fact, analyzing 
these nonlinear functionals very carefully, the explicit form of J-\^ implies that -1*1,2 are 
also C k ° mapping defined on IA. Rescale f = St, the explicit form of J-3 also implies \u% 
is a C k ° mapping in e, 5, r, and c. Substituting these observations into J^, it is straight 
forward to compute that h is a C k ° function with the expansion 

h(e,S,s,c) = c+ — + hi(e, S, s,c), hi = 0(S + e). 

47T 

Applying the Implicit Function again, we find that for small (e, 5, s), there exists a solution 
c = — j- + 0(e + 5) which is C k ° in e, 5, and s. Substituting this back to J-3 again, we also 
obtain u 3 = 0(S 2 + Se). □ 



5. Global bifurcation with a point vortex 

5.1. Introduction and notation. Consider the steady point vortex problem (PtV) in 
the periodic setting (Per) which corresponds to equations (2.5), (2.8), (2.7), and (2.12). 
Unlike in deriving small amplitude solutions in Section 3, we can not work with the e- 
scaled unknowns for the global bifurcation problem and it is more natural to work with 
the original unknowns (77, "0, c) , where the graph of T] defines the air- water interface, ip is 
the trace on the surface of the harmonic part of the stream function, and c is the wave 
speed. We have a point vortex situated at the origin with vortex strength e. The abstract 
problem is to find (e, rj, tp, c) in the zero-set of the operator 

where 

X := H^LS 1 ) x H^LS 1 ) x R, = {(r/,V,c) G X : 77(0) > -1}, 
Y := H^~ 2 (LS 1 ) x H^iLS 1 ) x R, k > 3 
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and 



(5.1) 




1 



T 3 (e; V ,^,c) :=c+(d X2 i; n )(0) + ^- 



k 2 7T 2 L 2 + l 



k=- 



oo 



Note that we applied the operator — d Xl to (2.8) to define T2, which is possible since we 
work on functions with zero mean on S 1 . The advantage of this procedure is to make T>^Ti 
a second order positive operator like D^T\. 

With our choice of spaces, J 7 is a C°° nonlinear mapping as long as the free surface 
does not contact the point vortex (in fact, T is analytic (cf. [7, Lemma 2.2]). Up to this 
point, it has been shown that there exists a curve of solutions ^i oc parameterized by the 
vorticity strength e bifurcating from the trivial solution e = 0. In this section we endeavor 
to continue ^\ oc globally. With that in mind, define Z C R x O to be the zero set of J- and 
let ^ be the connected component of Z containing the local curve ^ioc- Our main theorem 
is the following: 

Theorem 5.1 (Global bifurcation). One of the following alternatives must hold: 

(i) is unbounded inRx X ; 

(ii) there exists a nontrivial irrotational (i.e. e = 0) traveling wave solution in O where 
the velocity at the original is equal to the wave speed c; or 

(iii) along some sequence in {(e n ,ri n ,ip n ,c n )} C "io , we have r/ n (0) — > —1. 

Clearly Theorem 2.2 is an immediate consequence of the above statement, and so we 
need only concern ourselves with the proof of Theorem 5.1. The tool we elect to use is the 
global implicit function theorem of Kielhofer, which we now paraphrase. 

Definition 5.1. Let X and Y be (real) Banach spaces such that X is continuously em- 
bedded in Y. 

(a) A linear map A € C(X, Y) is said to be admissible provided that 

(i) A is a Fredholm operator of index 0, 

(ii) A : Y — > Y with domain of definition D{A) = X is closed, and 

(iii) the spectrum of A in a strip (—00, a) x (— ie, ie) C C for some a > consists of 



(b) An operator F = F(X,x) G C 2 (R x 0,Y), where O C X is open, is said to be 
admissible provided that 
(i) D x .F(Ao, £0) is admissible in the sense of (a) for each (Xq,xq) GlxO, and 



finitely many eigenvalues of finite algebraic multiplicity. Moreover, their total 
number is stable under small perturbations in the class of C(X, Y). 



WATER WAVES WITH COMPACT VORTICITY 



2V 



(ii) F is a proper map: for every compact U C Y, and K C X closed and bounded, 
K n F~ l (U) is compact in R x I. 

When F is admissible in the sense above, a degree theoretic global continuation argument 
(cf. [19]) can be applied to the local curve furnished by the implicit function theorem. 
Intuitively, the idea is that the principal part of the operator has a compact inverse, allowing 
a generalization of the classical Leray-Schauder degree. The result is a "global implicit 
function theorem." 

Theorem 5.2. ([14, Theorem II. 6.1]) Consider an admissible operator F = F(X,x) : 
K x O — > Y . Suppose that there exists a solution 

F(X ,x ) = 

at which D x F(Xq,xo) is an isomorphism. Let Z denote the zero-set of F and let c € be the 
connected component of Z that contains the local solution curve through (Ao,xo) furnished 
by the implicit function theorem. Then one of the following alternatives must hold: 

(i) <*f = {(A , x )} Uff+Uff-, where <tf+ ntf' = 0, and ^ + , <g- are each unbounded 
in M x X or with zero distance to M. x O c ; or 

(ii) \ {(Ao,xo)} is a connected set. 

We will apply Theorem 5.2 to obtain Theorem 5.1 so long as we are able to confirm that 
the operator J- defined by (5.1) meets the admissibility criteria of Definition 5.1. This task 
we address in the next section. 

5.2. Admissibility of T . The regularity of T is only guaranteed so long as the point 
vortex does not contact the free surface, and hence working near the boundary of O may 
pose some technical inconvenience. Instead, fix 8 > 0, and put 

O s := {(e,rj,^,c) elxO: dist(0, dn(rj)) > 5}, Z 5 :=Zr\0 5 . 

Let ^fg be the connected component of Z$ containing ^loc- All of our work in this section 
will be with Os and ^5. Indeed, it is precisely this restriction that is responsible for 
introducing the third alternative in Theorem 5.1; the first and second correspond to those 
in Theorem 5.2. 

Lemma 5.1 (Proper map). Let U C Y be compact and K C Ix Os be closed and bounded. 
Then Li n J-~ l {U) is compact inKxl. 

Proof. Let U and K be given as above. Put V := K n J 7_1 ([/) C K x O and let v n := 
{{(-niVni^niCn)} be a bounded sequence in V. By construction, 

— \fni 9m hn) ■ — •Fi^-ni Vni Ipni Cn) C U 

is bounded in Y. Compactness of U ensures that (modulo a subsequence) {u n } converges 
to some u in Y. In order to prove the lemma, we must use this information to construct 
a convergent subsequence of {v n } in V. S incG 6 n , c n G M, it is obvious that along b, 
subsequence e n — > e, c n — > c, for some e, c E M. 
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^From (5.1) and the definitions of u n , v n , we have 

fn = C n {G(Vn)4>n + e n (- V ' n , 1) • VG) + ^ {GiVn^n + 1) ■ VG)' 

- 2 ^ + ^ )2) W - VnGiVnHn + e n (l + (l]' n ) 2 )d Xl G) 

+ g(l + rj n )-b( ) + a 2 K(r] n ) 

n" 

(5.2) =: M(e n ,V„.,r?„,c n ) -a 2 - , " 3/2 , 

(1 + «) 2 ) 7 

where M is a bounded quantity in Y. Denoting 

M n = M(e n , tp n , rj n , c n ), M mn := M m - M n , K mn := K(r) m ) - n(r] n ), 
we see from above that 

fmn ■ — fm fn — M mn + OL K mn . 

Since {/ n } is convergent in H k ~ 2 (LS l ) by assumption, we know 

(5.3) d k ; 2 f mn ^0 in L^LS 1 ), 

Also, since {M n } is bounded in fl" fc_1 (LS 1 ) J the sequence {d Xl 2 M n } is bounded in i7 1 (L§ 1 ). 
Morrey's inequality implies {d k ~ 2 M n } is bounded in C°' 1 / 2 (LS 1 ), and hence it is precom- 
pact in C°' e (-LS 1 ), for any e £ [0, 1/2). This allows us to extract a convergent subsequence 
in C°' 1 / 4 (LS 1 ), say, which is in turn convergent in L 2 (LS 1 ). It follows that 

(5.4) fl£- 2 M mn ->-0 in L 2 (LS X ). 

From (5.3) and (5.4), we infer that d k ~ 2 K mn converges to in L 2 (LS X ). Meanwhile, 
because {r]' n } is bounded in H k ~ 1 (L§ 1 ), Up to a subsequence, therefore, (r)' n ) converges in 
^(LS 1 ), s < k — 1. Since k > 3, this observation permits the following estimate 

< l K mn|/ffc-2|(??m) 3 |/f fc - 2 + |K n |^fc-2 | (rj' m f - (l]' n ) 3 \ H k-a ->■ 

We have therefore proved that {r] n } has a convergent subsequence in H k (LS 1 ). 
Next we consider {ip n }- The definition of J- in (5.1) yields 

5n = -CnVn ~ 4>n ~ ^nd Xl ((1, ' VG) . 

We are given that {g n } is convergent in i? fc ~ 2 (L§ 1 ), and we have at this point demonstrated 
that {r] n } is precompact in H k (LS l ), as are {c n } and {e n } in R. The precompactness of 
VVi ^ follows immediately, which completes the lemma. □ 

We now compute the Frechet derivative of the operator J-(e; rj, tp, c) and proceed to 
prove the remaining admissibility properties. Recall G, VG, etc. are always evaluated at 
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X2 = 1 + rj(xi), we have at (e; rj, ip, c), 

Ji„C ={c + G(v)i> + e(-V, 1) • VG) (<S„(7?)C, <A> + 0) • VG 

+ eC(-V, 1) • Vd X2 G) +g( + a 2 K v (n)( - b v (e, rj, ij>, c)( 



+ T^iW- + e(l + (V) 2 )^G) 2 

(5.5) - t4s (^ - r/Om + e(l + MXiG) 

• R'£(^ - V(^WCV') + 2erj'C , d Xl G + e(l + (t/^X^G) 
^0 = (c + ^(r?)^ + e(-7/, 1) • VG)Q{r})4) - b v (e, 77, ip, c)4> 

" 7^2 - + + ( ) ?') 2 )9x 1 G) (</>' - r/^r?)^) 

•Fie = G{ri)i> + e(-rf, 1) • VG - b c (e, r?, V, c). 
Here the terms Db{- ••)(•) are only the average of the remaining terms in the expressions. 
72,,C = " ^ ((0, CO • VG + C(l, rf) ■ Vd X2 G) 

' > T A Al> T I' 

J-2xp4> = -<P , -T2C = ~f] 

(5.7) Fzn( = d X2 {n v (ri)CiP)(0), ^ = 8^(0), T 3c = l. 

From here on we shall denote by DJ 7 the Frechet gradient 

DF:= (D V T,D^,D C T). 
The next step is to show that DT is Fredholm. 

Lemma 5.2 (Fredholm map). At each (e,rj, V>,c) E Os, DJ r (e;rj,ip,c) is a Fredholm oper- 
ator of index 0. 

Proof. Fix (e, 77, ip, c) G Os- We seek to establish an estimate of the form: for all (£, tfi, d) 6 
X, 

(5-8) |(C,<Mlx < \D x T(e;n,M(<:,M\Y + \{CA,d)\ z , 

where Z is a space such that X is compactly embedded into Z, and the constant may 
depend on (e,rj,ip,c). An immediate consequence of (5.8) is that T is semi-Fredholm. By 
the connectedness of Os, and the fact that at the point of bifurcation DT is an isomorphism, 
we are then able to deduce that DF(e; 77, ip, c) is Fredholm of index 0. 

Estimate (5.8) is merely a consequence of elliptic regularity. To see this, up to constants 
to make the mean zero, rewrite J-i v equation of (5.5) as 

Ji„(e; r?, i\>, c)C = MiC + M 2 C' + Ma^faX, ^) + a%fa)C 
^(e; 77, 7/,, c )0 = M 4 0' + M 5 g(77)(/) 
J 7 i c {e;rj,ip,c)d = dM 6 . 
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Here the multipliers Mj belong to H k ~ 1 (LS 1 ) depending on (r],ip,c), and so the highest 
order term is that involving the derivative of the curvature operator k v : 

In light of the estimates in Lemma A.l, we may treat this as an elliptic problem: 

-C" = -^(v'fDTtiC&d) +.Mi(C,<M 

where Mi : Z ->• i? fc ~ 2 (LS 1 ) is a bounded operator and Z := H^iLS 1 ) x fl* -1 ^ 1 ) x R 
(which is indeed compactly embedded in X.). Prom this we have 

ICI//fc(L§i) i$ ICIl2(LS1) + I^X-^l(C<M)lHfc-2(L§i) + l-^i(C,^,d)|//fc-2( LS i) 

(5.9) < iDxTxiC^d^H^^ + KC^d)^. 
Next we turn to control of (f>. From (5.6) we have 

-<t>" = -^2r,(e; V, V>, c)C + T2^(e\ rj, tp, c)<f> + J" 2c (e; r?, ^, c)d 

+ c(" + ed xi ((0, C') • VG - eC(l, rf) ■ Vd X2 G) + rf'd 
=■ ^2r/(e; f?, V>, c)C + ^(e; ?7, c# + T 2c {e; n, ip, c)d + M 2 (C, d). 

It implies 

l^lff*(iSi) < l- D A-7 7 2(r?,V',c)(C,0,d)| Hfc -2 (L gi ) + |A^2(C,d)lH*-2(LSi)- 

Recalling the definition of M 2 , we have in fact that 

\M 2 (C,d)\ H k-2 {LS i ) < ICIh^ls 1 ) + 14 

and thus (5.9) implies 

(5.10) Mh^LS 1 ) ^ 1-^X^1 (Cj 0) c ^)l_H' fc - 2 (L§ 1 ) + \DxJ 7 2(C,<P,d)\ H k~2^ L gi^ ) + \(CA,d)\z. 
Lastly, in view of (5.7), 

d = D X T 3 (e; n, if>, c) (C, <P, d) - 8 X2 (H v ( V )(, V>> (0) - d X2 <foi (0) , 

whence, 

\d\ < \D x T 3 {e + 101 

This, together with (5.9)-(5.10), furnishes (5.8) and thereby the lemma. □ 

Lemma 5.3 (Spectral properties). The following statements hold. 

(i) For all 5 > and M, 9q > 0, there exist constants c±, 02,03 > such that, for all 
(e,ri,ijj,c) £ 0,5 wtt |e| + \ (r],ip,c)\x < M, and for all ((,<j),d) £ X, fi £ C with 
I arg/i| > | A* | > c\, k > I > 2, 

ICl/^(L§i) + l0ltf<(LSi) ^ c 2^l'" 2/2 (l^l(^)^l^-2(L§i) + |^2(/*)^U«-a(LSi) + \ d \) 

(5.11) , 

- C3^2\n\*(\(\ H l-j( L §i) + I^Ih'-^LS 1 ))- 
3=0 
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where 

I 1 

A := DxJ(e;i],)/!,c), A(/z) := A - (j, \ 10 

\ 1 

6ot/i of these being viewed as unbounded linear operators from Y to Y with domain 
of definition X and A\2 are the £ and (j) components in Y. 

(ii) For all (e, rj, if), c) G t/ie spectrum o~(A) consists of finitely many eigenvalues each 
of finite algebraic multiplicity, and has no finite accumulation points. Moreover, 
there exists an open strip S = S(e;rj,ip,c) around the negative real axis (— 00, 0] in 
the complex plane such that S n cr(A) is finite. 

Proof (i) Fix M, 6 > 0, fj, G C, \fi\ > 1, and (e,rj,i/j,c) elxlas above. Let £ = : 
R — )• R be a smooth cutoff function such that 

supp^C [-2,2], £ = Ion [-1,1], £>0, 

and put 

vf- = x ) := e -^l 1/2 'e(t)V;(x), ^ := (C, <t>, d). 
Let W:=Ix LS 1 . An useful observation for such type of functions is 

I 

(5.12) \e-^ 1/2t i{t)f{x)\w(u) < E H^I/lfl-W < |e-^l 1/2 ^(t)/(x)|^ (w) 

m=0 

Following Agmon [1], we augmented operator 



A:=A + 



...it) 



of 











Of 















where # := arg/i. Observe that, for j = 1,2, 

^ = e-^l 172 ^^)^-^ - |/i|e i9 e-^l 1/2t e(t)^ 

(5.13) + e^H^'VC*)^ - 2i\v\ l l 2 e ie e^' 2 ^' (t)^ 

= e-^ 1/2t d(t)M^ + e^-l'*! 172 )*^^) - 2i| M | 1 /2^(t))^., 

where ^1.2 is understood as £ or </>. 

We see that the highest-order terms in An are 

a 2 K v ( V ) + e l9 d 2 = -a 2 d Xl ((t/)" 3 ^) + e ld d 2 , 

which is a uniformly elliptic operator otvIA (as |0| > 6q and \ f]\H k {L&-) — M with k > 3). As 
in Lemma 5.2, we split A\ into the lower-order terms A4\, that are bounded in H k ~ 1 (U), 
and the higher-order terms above: 

(5.14) AiV = {e w d 2 - a 2 d xl {{r 1 ')-^d Xl ))^ 1 +Mi*, 



32 SHATAH, WALSH, AND ZENG 

where the linear operator Ai\ satisfies 

\Mi(C,<t>,d)\ H i(L§i) ^ ICIifi+i(isi) + Mtf«+i( iS i) + 14 k — l> l>0. 
Using (5.12), (5.13), and the elliptic estimates from (5.14) yields, for k > I > 2 

£ + |e" iH1/2i e(t)^i(/^)^|^-2(W) 

+ |e-^l 1/2 '(e"(t) - 2i|/i| 1/2 e / W)Cl^- 2(W) + l-Mi*!*-^ 

+ E H (i+1)/2 ici^-- w + E H i/a iA^iflW(z*) 

j=o i=0 

i-2 

(5.15) < \4 i -w 2 \Ai{m + \^\Q\L\L&)+Y.\v\ j,2 \^\HW(L^y 

3=0 

for |/i| sufficiently large. 

We now repeat this process with the second component of the linearized operator. The 
augmented problem becomes 

A 2 * = (~d 2 Xl + e ie d?)* 2 + M 2 ^> 

where the operator M. 2 satisfies 

\M 2 (C, <t>, d)\ H i m) < |C|^+ 2 (lsi) + \d\, k-2>l>0. 

Similarly the ellipticity of — d Xl + e l6 df and (5.13) imply, for k > I > 2, 

\*2\ Hm < M'- 2/ V 2 (/i)V^- w) + E \^ j+1)l Mm-^ {L ^) 

3=0 

1-2 

(5-16) +El^ /2 (KI^(L Sl ) + MI). 

3=0 

First multiplying a large constant to (5.15), adding it to (5.16), using (5.12), and then 
taking \fj>\ >> 1, we obtain, for k > I > 2, 

\(\hi(lsi) + \4>\h'(l§i) ^ \v\ 1 ~ 2/2 {\Ai(h)4>\ H '-2(lsi) + \M{ij)$\h>-i(l&) + \d\) 

i-i 

~ CE (KlH'-i(LSi) + |0U'-J(LS1))- 

i=o 

(ii) Note that A is a closed operator and an argument along the lines of Lemma 5.2 
shows that A(fj,) is Fredholm. Since A = A(0) has index 0, invariance of the index tells 
us that the same must be true for all [i. To establish that /x is in the resolvent set of A, 
therefore, it enough to check that A{jj) is injective. In fact, estimate (5.11) implies that 
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A{fj) is injective when \i satisfies the hypotheses of part (i). To see this, note that if (£, (j>, d) 
is in the null space of A(fi), where \x is such that (5.11) holds, then for k > / > 2, 

z-i 

(5.17) |Cltf«(LSi) + Mff'(LSi) < C2\v\ l ~ 2/2 \d\ - £3*52 \v\H\(\h<-3 (LSI) + l0|ff»-i(LSi))- 

3=0 

From the definitions, we have 

o = A 3 (/i)(c,<M = <M^)C,V>>(0) + 3^(0) + (l - fi)d 

which implies 

\d\ < i-rGClfl^S 1 ) + \^>\h^))- 
I A* I 

Substituting this into (5.17) and we obtain, for k > I > 2, 

l-l 

ICItf'(LSi) + Mff'(£Si) < ^l/ i l~(ICli/2(Si) + l0lH2(§i))-C 3 ^|At|^(|CI^^(L§i) + l^l^^(L§i))- 

j=0 

Taking I = 2 we obtain that <^>, and d has to vanish when \y\ >> 1 and hence the null 
space of A(fi) is trivial. We conclude that 

{z : | argz| > 9 , \z\ > c 2 } C cr(A) c . 

Furthermore, the compactness of the embedding ofXccF means that A has a compact 
resolvent, so that the spectrum is discrete, has finite multiplicity, and no finite accumulation 
points. These facts imply the second statement of lemma. □ 

5.3. Proof of main result. 

Proof of Theorem 5.1. Suppose there exists Sq > such that ^ C Os - Then Lemma 5.1, 
Lemma 5.2, and Lemma 5.3 show that J- is admissible in the sense of Definition 5.1 and 
thus Theorem 5.2 applies for all 5 G (0, <5o). Taking 5 = |5o, we see that in alternative (i) of 
Theorem 5.2, the solution curve can not approach the boundary of R x Os it it has to stay 
in M x Og , therefore it implies alternative (i) of Theorem 5.1. If alternative (ii) of Theorem 
5.2 happens, we claim it implies (ii) of Theorem 5.1. If otherwise happens, namely, there 
does not exists nontrivial zero point of F(e = 0, •), then ^ \ {(0, 0, 0, 0)} = ^ + U where 
<jp ± ■= <g n (M ± x O). Clearly C± is open iaV\ {(0, 0, 0, 0)} and thus <T \ {(0, 0, 0, 0)} is 
not connested. 

However, if there exists a sequence of 5 n > with 5 n — > for which ^ is not a subset 
of Os„ for any n, then there must be a sequence {(e n , r] n , <p n , c n )} C "5f and {x n } E E such 
that 

\(x n ,l + i] n (x n ))\ < S n . 

Clearly x n — > 0. If we assume that alternative (i) does not occur, then sup n \ r]' n \L°° ='■ M < 
oo. A simple calculation shows 

1(0,1 + ^(0))! < {\xn\ 2 + Hx n ) - nm 2 ) 112 + 5 n . 

< \x n \{M)+5 n . 
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It follows that rj n (0) — > — 1 as n — > oo, which implies the third alternative (iii) holds. □ 

Appendix A. Regularity lemmas 

Here we recall some basic facts about the regularity of the Dirichlet-to-Neumann and 
harmonic extension operators. Let us begin by defining the spaces in which we shall obtain 
estimates. For any domain !!cK 2 with H s boundary, s > 3/2, let 

1/1^.(0) := inf{|F|^ (R2) : F G H S (R 2 ), F\ a = /}. 

The homogeneous norm | • IjWq) is defined analogously. The spaces H S (S}) and H S (S}) are 

then defined as the quotient of H S (M. 2 ) and H S (M. 2 ), respectively with equivalence classes 
determined by the above norm. Likewise, we denote 

: = {/ G HS (V) ■ [ fds = 0}, H°(n) : = {/ G £T« (fi) : / is even in Xl } 

if f2 is bounded or symmetric in x±, respectively. The spaces H^LS 1 ), iJ|(R) are defined 
in the natural way. Using appropriately chosen local coordinates, it is possible to define 
analogous spaces for functions defined on dQ (cf. [20, §6]). 
Consider a domain £1 given by 

n = {x : x 2 < 1 + »7(si)}, Tj G H S (R) or ^(LS 1 ). 

We define the harmonic extension map T~L(r]) by 

n(rf)tp := ifu, 

where 

Atp H = in n, <PH\dQ = ^P- 

Notice that this relies on the existence of a well-defined trace operator in £(H s+1 / 2 (£l), H s (d£l)); 
the reader is again directed to [20] for a proof of this fact. The Dirichlet-to-Neumann op- 
erator J\f(rj) is then defined by 

M(7])ip ■= N ■ {Vip n )\ dn . 

The non-normalized Dirichlet-to-Neumann operator G{rj) is also often used. For any tp 
defined on R, let -0(^1,^(^1)) = defined on dQ. Then 

ds 

(g(7 ] )iP)(x 1 ) = ( V , )(x 1 )(M( V )^)(x 1 ,7 ] (x 1 )) = (M(t ] )iP)(x 1 , V (x 1 )) — 

where {rf) = + (v') 2 - is defined similarly for rj defined on LS 1 . The next lemma 
lists the regularity properties of these operators that we require. 

Lemma A.l. For ? ? G H s ° , s > 3/2, let H s := H S (R) or H^LS 1 ). Then the following 
statements are true. 

(a) The harmonic extension map U{rj) G C(H s {dQ), H s+1 / 2 (n))nC(H s (dn), iF + ^(ft)), 
s G (0, sq]. Moreover, its norm is bounded uniform in rj in a bounded set in H s ° . 
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(b) N{rj) G £{H s (dQ),H s ~ l (dQ)) and g(rj) G £(H s ,H s ~ l ), s G [1 - s ,s ], are self- 
adjoint and invertible. Moreover, their norm and the norm of their inverses are 
bounded uniform in rj in a bounded set in H s ° . 

(c) G{rj) is G°° in their dependence on rj. For example, the Frechet derivative of G{rj) 
with respect to r\ has the estimate: For any £ G H Sl , si(§, so] and s G [1 — si, Si], 

l(^(j]K,^)|^-i < ICIjr-iMtf's 

where the constant is uniformly controlled by \rj\ii s o- 

The proof of these properties is either straightforward or standard. In particular, the 
reader is directed to [22] and [20, §6]. 

Appendix B. Auxiliary lemmas for the vortex patch 

In this section, we prove some technical lemmas establishing the stated smoothness of 
the nonlocal operators used in the proof of Theorem 2.3 in section 4. 

Lemma B.l. Suppose that f3 G X s for s > 3/2, where X s is as defined in (4.10). Then 
r G H s+l / 2 (Bi), F G H s+3 / 2 (Bx) and the mappings 

(T(/3))(z)-z£C( y X s ,H s+1 ^(B 1 )), and (3 G X s h-» (F,a) G H s+3 ^ 2 (Bi) x R, 

is of class C ko+1 . 

Proof. Let j3 G X s be given and take 

r = T(/3) :=%05 + cosfl)+iC^(/3 + cos#), 

where T~L denotes the harmonic extension on B\, and C the harmonic conjugate on B\. It 
follows from the boundedness of these operators that T G H s+1 / 2 (Bi); it is easy to show 
additionally that T satisfies the assumptions in (4.5). 

Next consider F. Recall that the way in which it is defined is as follows: By the 
assumptions 7, we may solve the following problem uniquely 

AG = \d z T\ 2 ^(G) in B\, G G Hl{B x ). 

For 7 G C N ,N > s+k +l&nds> §, the mapping (P,G) i-> \d z T\ 2 -f(G) is C ko+1 from X s x 

3 1 

H2 +s (Bi) to H s ~z(Bi). The nondegeneracy assumption in (4.4), the Implicit Function 
Theorem, and the standard regularity theory imply the local existence and uniqueness of 
G G i? s+3/2 (,Bi) and that it is C h)+l in /3. By the positivity of 7, we know that AG > 0, 
and hence we may define 

a := [ AGdx, F = -G. 
JBi a 

Thus (F,a) solve (4.11). From this procedure, the lemma follows immediately. □ 

We remark that this implies that D and Do are each H s domains, as they are the image 
of B\ under T and 5T, respectively. 
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In the following, we consider the function VGlii, 1 + erj(xi)) . Given /3 G X s , it induces 
T and F which in turn lead to uj and uj as given in (4.2) and (4.12). Finally G is defined 
in (2.1). 

Lemma B.2. Let s > 3/2 and < e,S <C 1, the mapping (P,rf) \-> VG(-,1 + erj(-)) is 
C k +i j rom my b oun d e d set in X s x #| +1 (R) to H° +1 (R) x H s + l (WL). 

Proof. The vorticity is determined by (rj, f3) according to 

6 X 

where, recall, uj = A/, and D = 5T(B\). From Lemma B.l, we know that uj G H s ~ l l 2 {D) 
and C ko+1 in /3 hence, u G L 2 {R 2 ) n L°°(M 2 ) and is C ko+1 in f3. 

By construction, the vortex patch has approximate radius 5 and sits at the origin, while 
the free surface is a perturbation of the line {x2 = 1}, and the phantom point vortex in 
the air sits at (0, 2). For e and 5 sufficiently small, the strip 5:=Kx (1/2, 3/2) intersects 
neither D nor the point vortex but contains the entire free surface. One may verity through 
direct computation that, for any x',x" ^ S, 

log | • -x'\ - log | • -x"\ e H k {S) 

for any k > and is harmonic. Integrating this linear property, it is clear that the function 
mapping 

1 f oj(x'} 1 

ie54 G(x) = — / log|x - x'\ dx' log \x - 2e 2 | G R 

27r J^2 e 2ir 

is harmonic and belongs to H k (S) for any k > 0. Therefore, VG satisfies the same property 
and thus its trace VG(xi, 1 + erj(x\)), for 77 G H* +1 is in i/ s+1 (R). The symmetry in x\ 
and the C k ° +l smoothness of VG(-, 1 + £??(•)) follow from the symmetry of uj and rj and 
the smoothness of F in (3. □ 

Lemma B.3. Let s > 3/2, and define H by (4.13). Then H is of class C ko+1 (X s ; tf^S 1 )). 

Proof. Recall that 

ffCS) = ±-af log |r(.) - r(z)| 7 (iF(z))|r'(z)| 2 dz 

2vr y Bl a 

= / io g |r(.)-r(z)|AF(z)cte. 

In Lemma B.l it was shown that T G F s+1 / 2 (,Bi), F G H s+3 / 2 (B 1 ) L 
In the following we write H(j3) in a different form. Since £3 = A/, 

— log I • I * uj — f is harmonic in Do • 
2tt 
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Denoting by jVd the Dirichlet-to-Neumann operator on Dq, since f\dD Q = 0, it follows 
merely from the definition that 

1 | 1 ~ 

M? [(^log| • \*u)\ dDo ] = Viv[— log | -| *£-/], 

where Vjv[ - ] denotes the inner product of the outward unit normal on 8Dq with the trace 
of V[-] on dD Q . 

Similarly, as Zj is supported in Dq, 

— log | • | * uj log | • | G H 1 {Dq) is harmonic in Dq. 

Then, writing A/j)g for the Dirichlet-to-Neumann operator relative to Dq, we have 

log | • | * 2 - — log | • \)\ 8Dq ] = -V N [— log | • | * 2 - — log | • |]. 

Note that Vat here still refers to the outward unit normal for Dq, hence the minus sign. 
In the above, we have used that fact that uj G L 2 (M?) implies log | • | * uj G H 2 (M. 2 ) and 
thus V(log | • | * ui)\q£) G H^(8Dq) is well-defined. Summing these two identities yields 
the following 

■Yd (— log | • | * Z)\dD + .Yog (log I • I * s - ^ log | • |)| 9Do = lo S I ' I _ v ^/- 
In order to relate this to H, let us introduce the operator N(j3) ■ H S (S 1 ) — > H S ~ 1 (S 1 ) by 

M{P)(g) :=[N Do {goT{py l )]oT. 

Much as in Lemma A.l, M(f3) G {.(^'(S 1 ), H S '~ 1 (S 1 )), for all s' G [1 — s,s], it is C°° in 
f3 with uniform bounds for f3 in any bounded set in X s . For the complimentary domain, 
Dq, we introduce the analogous operator A/" c (/3) G C(H S {S 1 ),^ defined by 

jV c ((J)(g) :=[N D c{goT{py l )\oT. 

Lastly, to pullback the normal derivative on Dq to B±, let 

V N = V N ((3) : H S (B 1 ) -+ H^iS 1 ), 

be the operator 

where Vtv is the outward unit normal derivative on Dq. Then, because HoT -1 = =- log | • |* 
u;, and f = F o T" 1 , with some shuffling of terms we have 

(M + m c ){h - J- log |r|) = Vat(^- log |r| - f) - #-L log |r|. 

Z7T Z7T Z7T 

Since the right side of the above equation belongs to // S_1 (S 1 ), it is clear that H G iif^S 1 ). 
Moreover, the smoothness of N in V and the smoothness of V (and thus that of N) and -F 
in /3 imply the C 0+1 smoothness of H vn (3. □ 

Putting together these lemmas, we arrive at the following. 
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Lemma B.4. Let s > 3/2, and define T = (J-\, F 2 , F3) as in (4.16). Then T is well- 
defined in the sense that T{e, 5; 77, 93, c, 0) E y, for (e, <5; rj, cp, c, 0) £ M 2 x X . Moreover, 
T G C ko+1 (R 2 xX-,y). 

Notation. 

• For any scalar quantity a, we set (a) := y/l + a 2 . 

• ip% is the harmonic extension of <p to ^(77). When we wish to make the dependence 
on 77 more explicit, we shall write H(7])(p. 

• G{ji) = { r n')N('n), N(rj) being the Dirichlet-to-Neumann operator in the fluid do- 
main 

tt = Q(rj) ■= {(xi,x 2 ) : x 2 < tj(x{)} 

with vanishing at x 2 = —00. For regularity properties of G(rj), see Lemma A.l. 

• For a region f] C M. 2 , we denote the outward unit normal on Oft as N, and the 
outward normal derivative by V^. 

• k is the curvature of the free surface; viewed as a differential operator on r\ it is 

• If X and Y are Banach spaces and F = F(X, x) : Rxl — > Y is Frechet differentiable 
with respect to x, we denote its Frechet derivative evaluated at (Ao, Xq) by ^(Ao, %o) 
or D x F(X ,xo). 

• For Banach spaces X, Y, C(X\Y) designates the space of bounded linear operators 
from X to Y. 

• If X, Y are Banach spaces and L is a linear map from X to 7, we write ^Y{L) for 
the null space of L, and ffl(L) for the range. 

• When Xi, X 2 , and Y are Banach spaces and G = G(f, g) : X\ x X 2 — > Y is Frechet 
differentiable with respect to /, then we let 

(Gf(fo,9o)f, 9) 

denote the Frechet derivative of G in the direction /, evaluated at (fo,9o), an d 
applied to g. 
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